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In  previous  papers  It  has  been  shown  how  the  functional 
equation  technique  of  dynamic  programming  could  be  applied  to 
derive  variational  relations  for  Green's  functions.  In  this 
paper  we  apply  the  method  to  the  variational  problem  yielding 
the  equation 

(pu')'  -f  (r(x)  Xq(x))u  -  v(x),  u(a)  -  u(l)  -  0. 

The  introduction  of  the  parameter  X  enables  us  to  study  the 
resolvent  operator,  and  thus  to  derive  variational  relations 
for  the  characteristic  values  and  tne  characteristic  functions 
of  the  associated  Sturm-Llouvllle  equation. 
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FUNCTIONAL  EQ’JA'I'IONS  IN  THE  THEORY 
OF  DYNAMIC  PROQRAilMIKO — X: 

RESOU.TNTS,  CHARACTERISTIC  FUNCTIONS  AND  VALUES 

Richard  Re  11 man 
Sherman  Lehiman 


1.  Introduction 

In  previous  papers,  we  l.ave  sh.own  how  t[.e  functional 
equation  technique  of  dynaml'^  prot^ramnlng  can  le  applied  to 
derive  variational  relatione  for  kemelo  and  Green's  functions. 
In  [IJ.  the  Green's  function  associated  with  the  second  order 

equation 

(1)  u"  q(x)u  •  0,  u(a)  -  u(l)  •  0, 

was  dlecuaaed,  while  In  ,  analogous  meti-ods  were  appllea  to 
partial  differential  operators  to  obtain  the  Hadamard  varia¬ 
tional  formula.  In  [p]  ,  the  FredliOlm  Integral  'equation  was 
treated  by  similar  means,  and  Jacobi  matrices  were  discussed 

in  [>]  . 

In  this  paper,  we  wlsn  to  present  "ome  * 

•xtenalons  of  these  results.  Intro  iuclng  th-.e  narameter  X, 
we  consider  the  general  equations 

(2)  (p(x)u')'  (r(x)  Xq{x))u  -  0,  u(a)  -  0,  u(l)  >  0, 

obtaining,  as  In  the  papers  cited  al  ove ,  variational  equations 
for  the  resolvent  operator  as  a  function  of  a.  Utilizing  the 
■•romorphlc  nature  of  the  operator  as  a  function  of  X,  we  are 

able  In  this  way  to  derive  variational  equations  for  the 
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characterl Stic  values  and  fi-nctlons. 

Corresponding  results  are  obtained  for  the  vectorwnatrlx 
system 

(?)  (P(t)x')'  +  XQ(t)x  -  0. 

In  both  cases,  certain  assumptions  have  to  be  made  con¬ 
cerning  the  coefficient  functions,  q(x)  and  Q(t),  In  order 
to  be  able  to  consider  the  differential  equations  as  the  Euler 
equations  of  associated  variational  problems.  Since,  however, 
we  know,  from  the  results  of  Miller  aind  Schlffer,  ,  con— 
c«»mlng  Green's  functions  for  general  linear  differential 
operators  of  order  n,  (where  quite  different  methods  are 
employed)  and  from  corresponding  results  for  Fredholm  operators 
and  Jacobi  matrices,  that  the  relations  obtained  hold  under  far 
weaker  assumptions,  the  Interesting  problem  arises  of  deriving 
the  more  general  results  by  variational  techniques.  In  this 
paper,  we  present  a  method  of  analytic  continuation  which 
reduces  the  case  of  continuous  q(x)  to  that  of  positive 
continuous  q(x),  and  the  case  of  continuous  symmetric  Q(t) 
to  positive  definite  Q(t). 

In  a  separate  paper,  ,  we  have  sketched  ex^’enslons  of 
this  technique  of  ainalytlc  continuation  which  enable  us  to 
treat  non— selfadjolnt  differential  operators  and  non— 8ymr>e trie 
matrices  by  means  of  variational  methods. 

Finally,  In  a  brief  section  at  the  end  of  the  paper,  we 
Indicate  how  similar  methods  may  be  applied  to  obtain  corres¬ 
ponding  results  for  Fredholm  kernels  end  the  associated 


sharacterlatlc  values  and  functions. 
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2.  Variational  Prol len 

Consider  the  \  oundary  value  proble:.' 

(1)  (p(x)u' )  '  ^  q(x)u  -  v( X ) ,  a  ^  X  ^  1, 

u(a)  •  0,  u(l)  -f  au'(l)  -  0. 

We  shall  assume  that  p,  q,  and  v  continuous  functions 

on  the  closed  Interval  [a,l].  If  tr.e  cerresj- ondln;*,  Sturm- 
Llouvllle  problem  with  q  replaced  ' y  Iq  and  v  by  0 
does  not  have  1  as  a  ct  aracterlstlc  vailue,  then  the  unique 
solution  of  (1)  can  be  represented  In  the  forr. 

(2)  u(x)  K(  x,y  ,a  )v(y  )dy  . 

a 

The  fiuictlon  K(x,y,a)  Is  ailed  t' e  Green's  fun':'tlon  for  tie 
boundary  value  nroLlem  (l).  We  wls;  to  study  tbe  ierendenc'-- 
of  K  upo.i  a  by  means  of  tiie  fun'^tloml  equation  r.et;  od  of 

dynamic  programming . 

We  Imbed  (l)  in  the  system 

(5)  (p(x)u’ )  '  +  q(x)u  -  v(y  )  ,  3  v  x  v  , 

u(a)  ■  c,  u(l)  -f-  au'(l)  -  C. 

A  solution  of  this  system,  can  expressed  In  tie  form 
u(x)  -  Uq(x  )  +  C0(  X ) 


(4) 
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Where  1®  ®  solution  of  (l)  and  hence  equal  to  the 

right  side  of  (2)  while  ^(x)  Is  a  solution  of  tlie  system 

(5)  (p(x)^')'  q(x)^  -  0,  ^(a)  -  1,  ^(1)  -f  a^'(l)  -  0. 

We  also  consider  a  variational  problem  associated  with  the 
system  (3)#  the  problem  of  maximizing  J(u,v)  over  all  u  for 
which  u(a)  ■  c  where 

(6)  J(u,v)  (q(x)u^  -  p(x)u*^  -  2\/v(x))dx 

a 

-  (ud))^. 


This  variational  problem  has  the  property  that  a  u  yielding 
the  tnaxlnnim  Is  a  solution  of  the  system  (3).  To  prove  this, 
we  replace  u  by  u  In  where  t)(x)  Is  a  function 

such  that  T](a)  «  0.  We  obtain 

J(u  *'r),v)  -  J(u,v)  +  ^  (q'-in  ~  pu'T)’  -  VT])dx 


a 


u(1)ti(1)( 


,2l  /;l/  2  |2\. 

^  I  c/  vqq  -  pt'  )dx  - 


J 


In  order  to  have  u  yield  the  maximum,  the  term  must 

vanish  for  all  functions  q  for  which  T)(a)  -  0.  Hence 

y  ’  ^  (qu  (pu* ) 


I 


v)T)(x)dx  -  q(l)(u(l)  *►  au'(l)) 
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Consequently,  a  u  yielding  the  rnaxlmum  must  satisfy  the 
boundary  condition  at  1  and  must  satisfy  the  differential 
equation  In  the  Interior  of  the  inten'al  Qa.lJ.  Hence,  a 
solution  of  the  variational  problem  provides  a  solution  of  the 
boundary  value  problem  (3). 

In  order  to  use  the  variational  approach  one  must  make  an 
asauaption  on  p  and  q  sufficient  to  guarantee  the  existence 
of  a  maxlaium.  It  Is  sufflclent,  for  example,  to  assume  that 
p(x)  is  positive  on  the  closed  Interval  [a,l]]  and  that  the 
smallest  characteristic  value  of  the  Sturm— Llouvl lie  prob  em 

( pu • )  ’  "0,  u(a)-0,  u(l)-»-ou'(l)-0 

Is  larger  than  1.  If  q(x)  Is  uniformly  positive, 
q(x)  >  d  >  0,  over  0  <  x  <  a,  this  condition  holds  if  a  Is 
sufficiently  small,  or  If  d  Is  sufficiently  large.  TlUs 
assuBiptlon  guarantees  the  existence  of  a  unique  maximum  of 
J(u,r).  In  §7,  we  shall  show  by  analytic  continuation  how  the 
results  we  prove  can  be  freed  of  this  restrictive  hypothesis. 

At  the  moment  we  wamt  q(x)  to  be  positive  so  that  we  can  easily 
locate  a  region  within  which  no  ch.aracterlstlc  values  occur. 

3.  Dynamic  Programnlng  Approach 

We  now  let 

f (a ,c  )  -  max  J(u, v ) 
u(a)-c 

and  derive  a  partial  differential  equation  for  f  by  means  of 
A  technique  of  dynamic  programming.  We  regard  u  as  describing 
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a  policy.  The  variable  c  describes  the  state  of  the  system 
at  a.  The  result  of  following  the  policy  u  for  a  time 
interval  [a, a  -f  ^  Is  to  transfom  c  into  a  new  initial 
state  u(a  £i)  for  the  interval  [a  -f  /i,lj  .  According  to  the 
principle  of  optimality,  an  optimal  policy  u  over  the  interval 
[a,lj  with  initial  state  c  must  have  the  property  that  it  la 
an  optimal  policy  over  the  subinterval  [a  -f  a,i]  starting  with 
Initial  state  u(a  ^).  Translated  into  a  formula,  the 
principle  of  optimality  yields  the  equation 

f(a,c)  -  max  (f((a  -f  A),  u(a  -r  t^)) 
u(a)«c 

+  /**^  (q(x)u(x)^  -  p(x)u'(x)^  -  2u(x)v(x))dx.. 

We  proceed  formally,  assuming  that  f(-/,c)  has  continuous 
partial  derivatives  and  that  the  maximizing  u  has  a  continuous 
derivative.  Then  as  A 0, 

f(a  +  A,  u(a  +  ^))  -  r(a,c)  ♦  .^(a,c) 

>  ^(a,c)u'(a)A  +  o(A.). 

v>  C 

Consequently,  as  A. — >  0,  we  have 

“  j^{a,c)u'(a)  q(a)c^ 

u(a)-c  r"' 

-  p(a)u‘(a)^  -  2cv(e)  +  o(l)|  . 


TYie  quadratic 
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-  p(a)u'^ 


takes  Its  maximum  value  at 


.  1  ^ 
■  ?pTrr?T  <“■=' 


Hence  we  obtain  the  partial  differential  equation 


(1) 


^  (a,c)  -  — - —  -  2av(a)  +  o^q(a) 

^  ^p(a) 


4.  Variation  of  the  Qreen'a  Function 

Let  u  be  the  function  whlc!‘  ru-ixlnlzrr  J(u,y^  for 
given  a  and  c.  Py  usln^  (?.?)  ard  (2.^),  we  'an  find  an 
equation  which  connects  f(a,c)  with  tie  Green's  function 

K(x,y,a).  We  have 


(1) 


f(a,c)  fqu^  —  pu'‘  -  ?uv)v]x  -  (u(l)h 


a 


cT 

a 


^  u (qu  +  (  pu  *  )  '  —  V )  lx  —  "  uvdx  —  u(  1 )  I  ( 1  '  u  '  ( 


+  u (  a  )  t  ( a  )u  '  ( a  ;  -  (  u  ( 1  )  ) 


/)  1 


uvdx  :  I  (  c  )  u  '  ( a  )  -  ^  ^ 
a 


ill  u(l)t^(})  +  au'(lM 


•c/ 

a 


' ^  ‘  cv  ; X  (  a  )  ( '■ )  0  '  (  •  ' 


TJalng  the  fact  that  u^^  satlsflr-r  (h.')  ana  0  Ratlr.flor  (j.'), 
wo  obtain  by  Integration  y  rartr 


y  ^  Jirvdx  -  0((puy'  -f  qu^)d>: 

a  a 

u^((d^')'  -♦-  q0)j>.  '»■  0(l)p(l)uyi) 
a 

-  C5(a)L  (a  )u^(a) 

-  0'(:)p(l)u^(l) 

4-  o' (a  )r.  (a  )Uj^(a)  . 

Sin -e  {?(<'•)  “  -‘nJ  u  and  0  satisfy  tf.e  sane 

,ono£?onr'ou.'  louni'ir’y  'onilllon  ':t  1,  v/e  '^on-luie  at 

(')  0V  i>;  -  -  i  („)«,!.(. 0  +  (»('-);  (1  !U||j(l)  -  e(’(l)p(  1  )u,(l  1) 

=  -  T  (  a  a)  . 

'T’liu^  ,  til'  second  and  lilni  tt  on  ti,--*  iant  line  of  (l)  arc 
equal.  It  turTir.  out  Ic  'e  r,ov>^  onvcnl'»nt  to  uae  tl.c  lntey;i'al 
e  y  I' per  a  Ion.  We  f  ave* 

(‘O  f(i^  ^  u^vd>  ^  pvdy  +  '‘r(H)0'(a) 

a 

"  “  ^K(x  ,y  )v(  X  )v{y  )dxdy  -  2g  ^  ^^vdx  *  c^p(a)0'(a). 

"i  a  a 

We  olserve  tl.at  we  can  express  0  in  terms  of  a  partial 


derivative  of  K.  By  (2),  we  have 
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0(y)v(v)  !y  *  -  •-■  -  r  {  I )  ^  '  ^;.  (a,y,'^Jv(v' 


a 


for  all  contlnuo'j  funrt^or.' 


(M 


^  y?: 


-  -  r(a  ^ 

>»  /. 


Comt  lnlnc:  ^xrrei’clon  '/ v’-^n  '  y  (  '  ' 


/  .  .  *  4- 

•*»  ^  - 


the  partial  dlffor'cntl'i!  '-•quit  Jon  (  .2'.  w'  o  t''ln  ut^on 
equating  temr  Indopen.lc'nt  ol'  .■  , 


/\/'^  ^  (,.  .v,a)v(:,''  (:;)  ; 


a  a 


a 


-  P  (  1  )  i 


p  (  X  )  V  ( :■  )  i  / 


'  .1 


(a) 


Now  If  w?  "equate  coeJTlclontr”  of  ;•/.  o  tain  t 

relation 


^5) 


yK 

cJa 


( X  , ) 


1 


p(  a  ) 


C'(x^p(  ^  •  (a'  .  ..a2  -  I 


This  fomal  equating  of  -o*  J  n*  ■ 

use  of  the  synimetry  of  M 
such  as  continuity  of 


ij,  t  *  5  p  :  •,;  J  t  ■  t  .■ 


'> :  •  r:.'  ♦  J  o: 


3K  f  ^  :  /  W  ' 

rrr  (  x  ,  j.’  ,  t  '  — - 3  (  ^  g  ( 

i(a> 


We  shall  prove  tl  Jr  in  t  , 

In  the  consideration  of  t’ (jun.i”'  v;  .u--  •  -'o  l*''r  ( f  .  ' 
we  have  not  cons  lie  red  thp  ,  oun  ,  -  ■  or.  iJtlan.  ■  (  c: 
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u'(l)  -  0.  (Condition  u(l)  -  C  ''orre.Tnonis  forrrially  to 

a  -  0;  the  condition  u*(l)  -  0  correa:  onds  formally  to 
a  -  OD  In  t,^*  condition  u(l)  au  *  ( 1 )  -  0.  If  In  the  varia¬ 
tional  pioblem  of  rruixlr.lzlnr  J(u,v)  wc  orr.it  t^  e  terr 

and  add  tne  conrtralnt  u(l)  -  0  or  u’(l)  ■  0,  respectively, 
the  development  proceeds  ac  (..efore,  e>.-ept  all  teiT-s  which 
Involve  a  vanls...  Tlie  final  ronultr,  obtained  are  tie  same. 


Juntlf Icatlon  of  Equating  Coefficients 
Tie  forma]  proce.lure  used  in  th.e  derivation  of  (4.3)  -an 
ie  Justified  iy  t(,e  fellow  Inc;  lerjna. 

Le  iTta  .  I<e t  F C  x  , y  )  ^  li  contlnuour  function  on  th.e 
region  i  <  x  ^  ,  d  y  <  1  an  1  su;  p03 e  t i . a t  1’’ (  x  , y  )  -  P (  y  , x  )  . 


v 


1  v(  V  )v(  y  ^  ixdy 


n 


for  all  -ont  Inuou.  fund  lon^  •/ ,  ’__y  ..n  •  t  Ion  P  (  x  ,  y  )  -  0  , 

Froor.  .  Ir  t,  we  S'  ow  tt.  ;*  1  v  r.l.hea  on  t  c  ulaponal  , 
l.e.  ,  ?(",'')=  .  To  (io  t  Ir,  we  ta.a  •  nequen 'o  of  contlnuour 

functions  v^  wlti'.  *-  rosltlv*  an.  ’•  .ndlno  to  0,  eaci.  of 
wi.lci  vanir.i  es  Identic  T.iy  cutrlle  t:,t  Iritorvil  ^ 

ind  Integra  te.M  to  1  ov-a  *■  portion  >f  '  ~  t  ,  C  +  w  i^i, 

Is  In  tie  lnte,’'Vtl  •  rasslnr  tc  1 e  llrlt,  we  find  tiiat 
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To  prove  tl.at  F(r,n)  “  0  wtion  T  ft  r)  we  use  Tunctlonc 
with  t  <  'C  -  /c  which  Inteci-ate  to  1  over  the 

portion  of  ,  r  In  t;.e  Intrrval  ‘ind  also 

Integrate  to  1  over  the  portion  '^h  "r;  —  e  ,  e]  in 

but  which  vanish  identically  outr^  le  ther^'  two  Intervals. 
Letting  ^  — >  0,  we  find 

t  p(n.s')  >  F(r, :)  +  F(r^,Ti)  -  c. 

since  F(C,C)  •  F(r),r|)  -  0  and  F(T,t^)  -  F(-^,^),  we  have 

-  0. 

6.  Change  of  Varleble 

In  this  section,  we  shall  nvihe  a  ''’  anco  of  varlal  le  w'  1  . 

leads  to  a  different  expression  for  ^(x)  an  1  hen -e  an  alter¬ 
nate  expression  for  the  variation  of  tr e  Green's  function.  We 
make  the  chcuige  of  varlal  le 

( 1  -t  a  -  X )  ^ 
u  ■  CJ - -  w 

( 1  a  -  a  ) 

so  that  w(a)  "0  amd  w(l)  aw '  ( 1  1  -  C  wl  cn  ’.j  r-atln'le;; 
the  boundary  conditions  u(a)  -  ■,  u(])  -f  au '  (  1  )  -  0,  (This 
transformation  Is  valid  also  In  -a.se  a  -  C;  for  tro  oniltlon 

U’(l)  •  0 ,  1 .  e .  ,  a  -  OD  ,  we  St*  t  u  »  '  w . 

TV.en  by  (2.^  ) , 


p-i33a 


J(u,v)  -  J(w,v)  2c 


'?!  I  1  a  -  X 


q  ( X  )  w  (  X  )  d  X 


1 1  -f  a  —  a ; 
+  2  /■!  d>:  -  2= 


v(x)  dx 


]  a  -  a 


(^1  -*■  a  -  a) 


^  r ( 1 )^( 1 ) 

1  -f  a  -  a 


c  ■  < 


;  I  1  -♦•  a  — 


V 


(.1  •♦•  a  -  a; 


—  q  (  X  )  ;1  X 


1  — eJjlI - ,  dx 

( 1  4-  a  -  a  )*■ 


i  (^)° 


( 1  a  -  a  ) 


'^r  inn'cr-.ln?;  t:  ^  se-oni  lnt*ar;:'a]  by  lnt'?t^"3Mon  \  y  parts,  wo 
o^  t  a  1  n 


f 


J  (  U  ,  V  )  -  J  W  ,  V  - 


I 


1  a  -  X 


) 


-  q  ( >  ) 


-«■  a 


-  '1 

r  '  ( .  )  I  i 
-f  a  -  a 


+  F(a,c,7) 


■lore  F  In  In  lo;  on  ^ont  ol  u-  an  i  ii.  F'^n '< 

r 

!■(  i,  )  .  .'(,;,v)  -  '^.'X  ■- 

u(-  '»  K(a).<"  , 


1  -f  a  -  X 


-  ^ 


I  q  (  X  )  -  r  '  ( ) 


'1  ■*-  a  -  0  J 


+  y(i,  .v). 


■  •■.'i  X  l  ■  1  ;•  1  lu'  'V  1  '  '  '.’••n 

/I 


(  X  ,  y  ,  "  )  vv  ( •/ )  - 


.  -r  g  -  ,, 
_(  1  n  —  '\j 


-a 


q(y^  - 


a  —  a  I 


iy- 


n  ^  0 
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U  K(x,y  ,a  )v(y  )(jy  + 

a 


:  +  a  - 


+  a  - 


K(x.y,a)  I----  q(y) 

a  11  -f  a  -  a; 


f'  *  (  y  ' 

_» _ ijLx. 


+  a  -  .1 


dy 


Thu8  ,  by  ( 2  JO  , 


^(x)  - 


1  ♦  g  —  X 
1  -f  a  -  a 


J 


K  (  >;  ,  y  ,  ;i  ) 


12  a  -  a ; 


_  !■'(}) 


1  -f  a  -  a 


Substituting  thlG  expression  Into  ( •• .  0  '■•e  rind  anoti.er  ror’^’-ula 
for  the  variation  of  tf.e  Green 'r  fun  t!  on, 


(1)  .  _1_ 

O’  a  p  ( a  ) 


1  ♦  g  -  X 

1  +  a  -  a 


o  1 


‘  K(x,v,a) 


1  ♦  Q  - 
1 2  u  - 


q(.v) 


i_LLi 


IVV 


.  +  a 


1  a  - 


a  - 


K  (  .  .■  ,  a  ) 


2 _ 

4  a 


--  q(>  ' 

'^J 


.  a  - 


j 


TTie  cases  u(l)-C  and  a;'*'  ■•anlj  ■  .nJ.  ei.  'y.o 

refults  are  what  would  he  oitcilneJ  ■‘♦vir  ■  a  -  f  on 

letting  a  — >  oo  res; ec t Ively  . 


7.  Analytic  Continuation 

In  order  to  extend  t^.e  foregolnp:  approacli  to  the  gener'a2 
cese ,  we  shall  employ  the  method  of  analytic  "'ontlnuat  Ion . 

Consider  the  equation 
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(1)  (pu-  )'  +  (e  >  q(x))u  -  ▼, 
u(a)  -  0,  u(l)  +  au'(l)  -  0, 

whtr*  B  la  a  poaltlva  quantity  ohoaan  larga  anough  ao  that 
2  -*■  q(x)  2  d  >  0  In  [x,l]. 

Pu2*thanBora4  wa  know  from  tha  Sturmlan  oonpax*laon  theoram 
that  z  oan  ba  ohoaan  la2*ga  anou^  ao  that  A  -  1  la  not  a 
chax*actariatio  valua  of  tha  SturiHLiouvilla  problam 

(2)  (pu'  )'  +  A(z  >  q(x))u  -  0, 
u(a)  «  0,  u(l)  au'(l)  «  0. 

It  followa  that  tha  inhOBOganaoua  equation 

(3)  (pu-  )'  (z  q(x))u  -  V, 

u(a)  -  0,  u(l)  +  au'(l)  -  0, 

will  than  hava  a  unique  aolution. 

The  Graen'd  function  aaaooiated  with  thia  problam  will  ba 
a  function  of  z.  If  we  oan  show  that  thia  function  is  a 
maroDorphio  function  of  z,  it  will  follow  that  tha  relations 
originally  obtainad  undar  tha  asaumption  that  z  is 
■uffioiantly  large  will  actually  ba  valid  for  all  z  diatinot 
from  a  aat  of  polaa. 

In  particular,  tha  ralationa  will  ba  valid  for  z  -  0, 
providad  that  ^  -  1  ia  not  a  charaotarlatic  value  of  tha 
equation 
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(4)  (pu-  ) '  q(x)u  =  V, 

a)  -  c,  u(l)  >  au'(l)  0. 

Oenerallj  speaking,  the  relations  will  be  valid  whenever  they 

■ake  sense. 

8.  Analytic  Character  of  Preen  s  Function 

Although  It  follows  from  well-knovm  results  in  the  theory 
of  linear  differential  equations  that  the  Green's  function  Is 
a  laeroBttrphlo  function  of  z,  we  shall  outline  a  proof  here 
for  the  sake  of  conpleteness. 

'fo  simplify  the  notation  let  us  write  o(u)  «  0  for  the 
tMundary  condition  u(l)  +  au'(l)  =  0.  Furthermore  we 
observe  that  the  following  considerations  also  cover  the 
boundary  condition  u'(l)  0,  so  that  we  could  also  take 

a(u)  -  u'(l). 

To  begin  with,  the  two  solutions  of 

(1)  (pu' ) '  +  (z  +  q(x) )u  -  0, 
detenslned  by  the  Initial  conditions 

(2)  u^(a)  -  1,  u^(a)  -  0, 


U2(a)  -  0,  u^(a)  -  1, 


are  entire  funotions  of  z  for  x  in 
fixed  X,  their  derivatives  u^  and 
funotions  of  z.  Hence,  In  particular, 


0,1  .  Also  for 
u^  are  entire 

0(02)  Is  an  entire 
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function  of  z  which  doct  not  vanish  Identically. 

Concerning  p(x)  and  q(x),  we  aur«  assuming,  as  before, 
that  p(x)  2  (3]^  >  0  In  [o,l  ,  that  It  has  an  Integrable 
derivative,  and  that  q(x)  la  Integrable  In  |^a,l  ,  and 
uniformly  bounded. 

The  general  solution  of 

(3)  (pu'  )'  +  (z  +  q(x))u  -  V 
can  be  written 

(4)  u  -  c,u.(x)  -*■  c^u  Jx) - v(y)dy. 

p(a)  a  L 

Imposing  the  boundary  conditions 
u(a)  0,  o(u)  ^  0, 

we  see  that  this  particular  solution  has  the  form 

u  -  ^  K(x,y  )v(y  )dy 


where  K(x,y),  th^  deelred  Green's  function,  is  given  by 


K(x,y ) 


r(a ) 


-  U2(x),  a  <  y  <  X, 


r 

I  io(u^)u^(x)  -  u_j^(x)o(a^) 

hTtL 


u^Cy),  a  <  X  <  y. 


We  see  that  for  any  x  and  y  In  a,ll,  K(x,y)  and  its 
partial  derivatives  with  respect  to  x,y  and  a  are  mere- 
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■orphle  funotiens  of  i.  iliac,  we  observe  that  thin  formula 
provides  a  proof  for  the  a/n&etry  of  the  Oreen's  function,  a 
property  which  we  have  used  in  the  previous  sections. 


9»  Alternate  Derivation  of  Expression  for  ^(x) 

The  explicit  representation  for  K(x,y)  given  above 
shows  that 


dK 


1 

x-a  p(a) 

It  followa  then  that  —  p(a)^(a,y,a) 
two- point  boundary  value  problem 


Is  the  solution  of  the 


(pu')'  qu  -  0,  u(a)  »  1,  o(u)  =  0. 


Hence  we  have 

¥(y)  -  -  P(a>^(a,y,a). 


The  extension  of  the  relation  In  (1)  is  the  keystone  of 


the  treatment  oi  Kl^lor  and  Schlffer  In  6|  ,  and  the  analogue 
of  this  relatlcr  for  partial  differential  operators  Is 
essential  in  thr  derivation  of  the  liadainard  variational  formula 
given  in  [4] . 


10,  Variation  of  Charaoterietlc  Values  and  Characteristic  Punctions 
Consider  the  Stum- ULouvl lie  problem 


(1)  (pu')'  +  (q(x)  Ar<x))u  -  v(x),  a  <  x  <  1, 

u(a)  -  0,  u(l )  ■♦'  au  ' (1 }  -  0. 


ItovlMd 


This  problen  i  aa  a  unique  solution,  \'heD  X  Is  not  a  charac¬ 
teristic  value,  which  Is  given  'v 


u(x) 


^  ( X  ,  y  ,  X  ,  a  )  V  ( y  ^  dy . 


e  function  called  the  resolvent,  lo  a  neromorphlc 

function  of  X  with  pcles  at  the  chai’acterlstlc  values 
X-  ,  X- ,  X  ,  ....  If  (x)  Is  t!.e  characteristic  function 

1  rz  ^  A 

inso'lated  with  t!.e  >^1  luractc rlstl value  X^^ ,  then  H  has 
the  re:  r’^rcntatlon 


(  X  j  '-.  ( y ) 

I  K 


:  ro:  lc’'i  (1)  Is  ol  tainod  f ror.  (f.l)  ly  replacing  q(x) 
•y  q(x^  Ar(;''>.  -'ror  (^.')  we  or  tain  the  equation 


( 


x\  ^  ( X , y , > , a 
1  a 


/  \  /  \  /  \ 
r(a)  —  (a,y,X,al  ^  (x,a,X,a). 
>  /  y  f 


"his  equation  Is  valid  for  ail  X  f  X.  wlth.out  any  ossur.ptlons 

t 

on  q(x)  an(i  r(:>:)  ex''e’  t  'onllni.lty  •  a.s  ^  eon  shown  ly 


analytic  continuation  In  ;V. 

Corr.l  Inlng  (2)  with  (3),  we  find 


f 


00  ( y  )  X, . 

r:*  X _ K  _ ^ 

r-1  “(X  -  X, -a 


(v) 


o'  a 


I 


J 


( 


00 


p(a) I  5 

VK-1 


'^w(»)'t,.(y )"( 1  qo  'ii,(x)V' (a)"! 


X  —  X, 


vX»l 


-  \  ; 
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Let  ting  X  — >  and  equating 

i' 

find 


■oelT  1 '' lents  of  ( 


(M 


JaT 


r(a)(V  (a))". 


I 

By  equating  coefflclentc  of  (x  -  \,  )  ‘  find 


(5) 


a\(x) 

“Xa~ 


p(a)'^  (a)  ^  i - 

^  V,  X  -  > 


We  can  also  derive  '^x:  rosr 'on~  :  o"  Mie  varL 


sums 


S„(x,y,a)  .  2 


00  (x)'*' (y) 

-  ,  (n  -  r,  3, 

k-1  K 


Since 


1 


1  X  X" 

T - ? ^ 


we  have 


00 

R(x,y,X,a)  -  -  2  ^  ^ 

n«l 


n-1 


Consequently , 


PS 

(x,y,a)  - 


-  p( 


n)  2 


•n  +  1 

i  <  J 

I?!-'  <n 


>s,  ( 


1 1 Ion  of  t  : 


\ 


y 
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11.  Matrix  Case 

Analogous  results  oan  be  obtained  for  a  Oreen's  function 
associated  with  the  vector-ouitrlx  system 

(1)  ^  ^  ^  ^  ^  Q(t)x(t)  -  y(t),  a  <  t  <  1, 

x(a)  =*  c,  x(l)  Hx(l)  -  0, 

where  P,  Q,  and  6  are  n  >*  n  matrices  and  x,  y,  and  o 
are  n-llnenslonal  vectors,  infe  shall  assume  that  ?,  Q,  and 

r 

y  are  continuous  on  the  closed  lntez*val  i^*lj  that 

P(t)  has  an  Inverse  for  each  t.  Also  we  shall  aosuioe  that 
P,  Q,  and  P(i)B  are  syimaetrlc  matrices  so  that  the  system 
Is  self-adJolnt . 

Por  notation,  we  denote  the  tra.  %  ose  of  a  matrix  A  by 
A*  and  we  use  (u,v)  to  denote  the  Inner  product  of  two 
vectors  u  and  v. 

Associated  with  this  boundary  value  problen  there  Is  the 
variational!  problem  of  maximizing  J(x,y)  wher^ 

J(x,y)  -  {PJi,x)  -  2(y,x)jdt 

-  (P(1)BX(1),X(1)) 

subject  to  the  condition  that  x(a)  -  c.  Let 

f (s,c )  max  J(x,y ) . 
x(a)-c 


As  In  the  scalar  case,  a  vector  x(t)  which  yields  the  maxi- 
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mua  otust  be  a  solution  of  the  boundary  value  problem  (1). 

Proceeding  as  In  one  can  use  the  principle  of 
optlsuillty  from  dynamic  prograamlng  to  derive  the  partial 
differential  equation 

(2)  -  |;(»,o)  ■  -  2(c,y(a))  MtKajo.o) 

where  If  .  . .  |C^  are  the  com^^onents  of  c  then 


We  wlah  to  use  this  partial  differential  equation  to 
study  the  Oreen's  function  for  the  system  (1).  The  Qreen  a 
function  Is  an  n  >  n  matrix  K(t,3)  for  which 

(3)  K(t,8)y(8)ds 

is  a  solution  of  (1)  with  *q(®)  ’  before  a  solution 

of  (1)  with  x(a)  «  c  can  be  written  as  a  sum 

x(t )  -  XQ(t)  +  ^(t )c 

where  ^(t)  is  the  n  n  nvatrlx  fujictlon  which  Is  the 


solution  of 
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l>{a)  --^  I.  -5(1;  ^  Bi<l)  ^  0. 

Sliice  x(i)  lE  the  maxi  ml  zing  vector  for  the  variational 
problem,  we  obtain 

(3)  f(a,G)  (y(t  ),x(t )  )dt  -  (  p(  1  )b(  1  )x(  1 )  ,x(  1 ) ) 


>  (P(a)x(a),x(a))  -  ( P( 1  )x( 1 ) ,x( 1 ) ) 

-  (y(t  ),x  (t )  )dt  (y(t  ),i(  t  )c  )dt 

^  \  a 

4  (r(a)xQ(a),c  )  4  (P(a)5(a)c,c  ). 


Also 


(-) 


(y  ( t ^  ( * '  ( t )  ,-^r(  t  )5(  t  )c  4  qJc  )dt 

a  a 


4  (r(l  )x^(l  ),5(l)c  )  -  (  p(a)xQ(a),i(a)c  ) 

-  (x^(l  ),P(1  )f(l  )c  )  4  (xQ(a),  P(a)i(t)c  ) 

-  -  (?(a)x^(a),c)  -  (F(l  )3Xq(1  )5(l)c  ) 

-  -  (p(a)xQ(a),c) 


because  ol  the  assumption  that  P(1)B  is  symmetric. 
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Hence  by  (3)  and  (5), 

f(a,c)  r,-.  (K(t,8)y(a),y(t )  )d8dt 

a  a 

-  2^"^^  (y(  t  },5(t  )c  )dt  ^  (r(a  Vf(a)c,c  ). 

a 

klLO  It  Is  possible  to  express  J  in  terms  of  a  partial 
derivative  of  K.  3y  (o), 

(c,i*  (  8  )y(  a  )  )dt  -  _  (c,F(a)^/'^^  -^(t,3,a)!  y(8)d8) 

a  a.  t=a 

for  all  continuous  vector  functions  y  and  all  vectors  c. 

Henoe, 

(7)  V('')  *  -  P(a>#(t,.,a)|  . 

^  I  t-a 

Next,  we  combine  the  above  expression  for  f(a,c)  with 
the  partial  differential  equation  (2)  a..d  equate  the  terms 
independent  of  c  to  obtain 

r,l  (^^y“.i.Qu^,iy(B  ),y(  L  ;  ^dsJt 

a  s 

a  a 


or 


Now  we  "equate  coefficients"  to  conclude  that 


(9) 
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ThlB  formal  argument  can  be  Juatlfled  by  the  fo- lowing  reault 
correapondlng  to  •:he  lemna  proved  In  ^8* 

Lemma .  Let  M(  t ,  a )  be  a  continuoua  n  ^  n  matrix 
function  for  and  auppose  that  M  has 

th.e  syr.unetry  property 

M^(t,3)  =  M(u,t). 


I£ 

(M(t,e)y(a),y(t  ))d3dt  -  0 

a  a 


for  all  cof'tlnuoua  vector  functloaa  y,  then  M(s,t)  ^  0. 

We  omit  a  proof  of  thla  lemma  becauae  one  can  be  oon- 
atructed  In  a  way  quite  similar  to  that  employed  In  It 

la  al3j  rot  difficult  to  verify  that  th'  mat’'lx  function  In 
(6)  haa  the  I'equlrel  ayrnmetry  pixjperty  teoause  the  system  la 
aelf-adjolnt . 

Combining  (  ,’ )  nnd  (y)  wo  obtain 


(10) 


yK 

9 


t ,  a ,  a  )  ^  t ,  8 ,  a  )  F(  a  )|j(  t ,  0 ,  a  j  1 

a^a  I  t  =>a 


The  problem  we  have  Just  l  j^eated  f  oes  not  certain  as  a 
special  case  the  problem  with  the  boundary  condition  x(l)  =  0. 
However,  o:\ly  a  small  change  Is  required  to  handle  this 
problem  by  the  sane  method.  In  the  definition  of  J(x,y)  one 
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oan  omit  th«  ttnn  —  (  F(  1  )B(  1  )x(  1 ) ,  x(  1 ) )  and  for  the  maxl- 
mlwitlon  problem  add  the  conetralnt  x(l)  =  0.  The  remainder 
of  the  argument  la  quite  aimilar  ai  d  the  final  result 
obtained,  equation  (10),  is  the  same. 

The  approach  by  way  of  a  change  of  variable  as  in  can 
also  be  followed  in  the  matrix  case.  Set 

x(t)  =  (I  >  (1  -  t)Bj(l  +  (1  -  a)B)~^c  +  w(t). 

•mie  transformation  has  the  propeirty  that  if  x(a)  c  then 
if(a)  -  0;  and  if  x(l)  +  Bx(l)  =  0,  then  w(l)  Bw(l)  0. 
Here  we  soist  also  maice  the  additional  assumption  that  the 
inverse  (I  (1  —  a)B)”^  exists.  Proceeding  as  in  ,  we 

obtain  the  equation 

/ 

J(t)  .  I  +  (1  -  t)B  K(t,»)|Q(8)(I  +  (1  -  8)D) 

a  ^ 

-^(P(t))Bjd 

•  ( i  ■»  (1  -  a)B)“^, 

which  allows  the  derivation  of  another  expreselon  for 
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12.  Integral  Equations 

Let  U8  now  Indicate  briefly  how  the  same  formallam  may  be 
applied  to  Integral  equations. 

The  equation 

(1)  u(x)  -  v(x)  4  k(x,y)u(y)dy 

a 

has  the  solution 

(2)  u(x)  -  v(x)  K(x,y,X)v(y)dy, 

a 

where  K(x,y,x)  Is  the  Fredholm  resolvent. 

If  k(x,y)  Is  positive  definite,  we  can  consider  (l)  to 
be  the  Euler  equation  corresponding  to  the  problem  of  minimi¬ 
zing  the  quadratic  functional 

(3)  Q(u)  -  k(x,y)u(x)u(y)dxdy 

a  a 

+  2 U(x)v(x)<l3t  u®(x)dx. 

a  ‘a 

Regarding  the  minimum  over  u  as  a  function  of  a  and 
functional  of  v,  we  can  proceed  very  much  as  above.  A 
derivation  of  the  variational  equation  for  Fredholm  resolvent 
Is  given  In  ,  Continuing  as  In  we  can  derive 

variational  equations  for  the  characteristic  values  and 
functions. 
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